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Abstract. We investigate bright and dark solitons with anomalous or normal 
dispersion and under transverse harmonic confinement. In matter waves, positive 
atomic mass implies anomalous dispersion (kinetic spreading) while negative mass 
gives normal dispersion (kinetic shrinking). We find that, contrary to the strictly one- 
dimensional case, the axial and transverse profiles of these solitons crucially depend 
on the strength of the nonlincarity and on their dispersive properties. In particular, 
we show that, like bright solitons with anomalous dispersion, also dark solitons with 
normal dispersion disappear at a critical axial density. Our predictions are useful for 
the study of atomic matter waves in Bosc-Einstein condensates and also for optical 
bullets in inhomogeneous Kerr media. 



PACS numbers: PACS Numbers: 03.75.Lm; 42.65.Tg 
1. Introduction 

Matter waves made of Bose condensed alkali-metal atoms [1] and optical waves in 
nonlinear Kerr media j2] are accurately described by a three-dimensional nonlinear 
Schrodinger equation (3D NLSE). In fact, the NLSE is a unifying theoretical tool of 
nonlinear optics and of the new field of research called nonlinear atom optics |H1 Hj. 
Many experiments have been devoted to the study of solitary waves in both matter 
waves EH 13 IHl M an d optical waves [TP] . Bright and dark solitons are usually 
analyzed in quasi-lD configurations, where analytical solutions of the ID cubic NLSE 
are available fUE]- O n the other hand, it has been shown that dimensional effects can 
be quite important: under transverse confinement atomic bright solitons with anomalous 
dispersion, i.e. kinetic spreading, become unstable if the transverse confinement is not 
sufficiently strong [T3]. 

In this paper we analyze bright and dark solitons with anomalous and normal 
dispersion (kinetic shrinking) taking into account the transverse width of the solitary 
wave. We introduce a non-polynomial Schrodinger equation, which extends that we 
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derived some years ago by including the case of the normal dispersion. In this 
way, we find new analytical solutions for 3D bright solitons under transverse harmonic 
confinement and show that, contrary to the strictly ID case, 3D black and gray solitons 
with normal dispersion exist only below a critical axial density. Our predictions can be 
tested not only with matter waves in a periodic optical potential along the axial direction 
[To] but also with optical bullets in a inhomogeneous graded-index Kerr medium fTol . 

2. Matter waves and optical pulses 

The 3D cubic nonlinear Schrodinger equation (3D CNLSE), which describes both matter 
waves in a Bose-Einstein condensate and optical pulses in a inhomogeneous nonlinear 
Kerr medium, is given by 



where g > and the parameters 8d = ±1 and 5jy = ±1, according as whether the axial 
kinetic term has anomalous (5p = 1) or normal (5d = —1) dispersion, and the nonlinear 
cubic term is self-focusing (5n = 1) or self-defocusing (5n = —1). In the context of 
nonlinear atom optics, the field ip{r,t) of Eq. (1) is the macroscopic wave function of 
the Bose-Einstein condensate. In Eq. (1) the time t is in units of 1/uj±, where uj± 
is the frequency of the transverse harmonic potential. The lengths are in units of the 



characteristic harmonic length a± = yK/(mu±), where m is the atomic mass. The 
mass m is the true atomic mass in the absence of an axial confinement (Sd = 1) [T3lll4| . 
but it is instead the modulus of a negative effective mass which takes into account the 
effect of periodic axial potential when atoms are around the upper band edge (5r> = —1) 
|15| ITo] ITTj. The strength g of the cubic nonlinearity is proportional to the modulus of 
the s-wave scattering length of the inter-atomic potential and its sign is given by the 
parameter 5n'- o~n = 1 corresponds to an effective interatomic attraction and Sn = — 1 
corresponds to an effective repulsion [TB^ ITU [TU] . 

In the context of nonlinear guided wave optics, the field ^(r, t) of Eq. (1) represents 
instead the envelope of the electric field oscillating at a fixed frequency. Both the cubic 
nonlinearity and the transverse harmonic confinement model the refrective index of an 
inhomogeneous Kerr medium Note that in describing an optical pulse in fibers with 
Eq. (1) the two variables t and z have an inverted meaning: the variable t is the axial 
coordinate of the electric field while the variable z is the time coordinate of the electric 
field. The scaling of these variables is discussed, for instance, in Ref. |17j . where g is 
set equal to l/(2ir) and the field ip(r, t) is normalized such that / \ip(r, i)| 2 ci 3 r represents 
the constant pulse energy. 

To study the 3D CNLSE time- dependent variational methods are often used 
[21 UH E21 ^] • The 3D CNLSE is the Euler-Lagrange equation obtained by minimizing 
the following action functional 
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Some years ago we proposed a variational ansatz of the field ip(r, t) for the Eqs. (1) 
and (2) with Sd = 1. The resulting effective non-polynomial Schrodinger equation is very 
accurate in reproducing the numerical results of the 3D CNLSE with Sd = 1 [T^| I19j. 
This effective equation has been used also to model the data of various experiments 
(2Hj. Here we use the same approach by considering also the case with Sd = —1- The 
variational ansatz is 

^M) = f(z,t) YTT-T.. ,, » ( r /M 



' x 2 + y 2 ' 



'^/ 2 a(z,t) 2 ^\2a(z,t) 2 J ' (3) 

where f{z,t) is the axial wave function and o~(z,t) is the width of the transverse 
Gaussian wave function. By inserting Eq. (3) into Eq. (2) and integrating along x 
and y, the resulting effective action functional depends on two fields: f(z, t) and o~(z, t). 
Neglecting the derivatives of o~(z,t) we find that the Euler-Lagrange equation of the 
axial wavefunction f(z,t) given by 



.df 

''m 



S_d&_ 

2 dz 2 



S 



+ 



a- 



f, 



(4) 



and the Euler-Lagrange equation of the transverse width o~(z, t) is instead given by 



a 



-S N g\f\ 2 . (5) 

Inserting this formula into Eq. (4), we obtain a nonpolynomial Schrodinger equation 
(NPSE). In the weak-coupling limit g\f\ 2 <C 1, where a ~ 1, the NPSE reduces to the 
familiar ID cubic nonlinear Schrodinger equation (ID CNLSE) 



.df 

''dt 



2 dz< 



-S N g\f\' 



(6) 



where the additive constant 1 has been omitted because it does not affect the dynamics. 

It is well known that the ID CNLSE admits solitonic solutions. In particular, 
for 5d ■ o~n = 1, i-e. for anomalous dispersion (6~d = 1) and self-focusing nonlinearity 
(5n = 1) but also for normal dispersion {5d = — 1) and self-defocusing nonlinearity 
(5jv = — 1), one finds the localized bright soliton solution 

g(z - vt) 



f(z,t) 



—- Sech 
4 



i5 D v(z-vt) e i(v 2 /(2S D )-fi)t 



(7) 



where Sech[x] is the hyperbolic secant, v is the arbitrary velocity of propagation of 
the shape-invariant bright soliton and the parameter /i = —Sng 2 /^ is obtained by the 
normalization condition of the field f(z,t) to one fH\ I12j. Instead, for 5e> ■ <5jv = — 1) 
i.e. for normal dispersion (Sp — — 1) and self-focusing nonlinearity (5^ = 1) but also 
for anomalous dispersion (5d = 1) and self-defocusing nonlinearity (5^ = —1), one finds 
the localized but not normalized dark soliton solution 



f{z,t) 




Tanh 



(z-vt)Jg<j) 



+i 
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where Tanh[x] is the hyperbolic tangent, (p^ > is the constant value of the field 



amplitude \f(z,t)\ at infinity 



±oo), and v is the velocity of propagation of the 
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dark soliton. The parameter /i = Sogcjf^ is obtained by the asymptotic behavior of the 
field f(z,t). For < \v\ < c s , where c s = Jg^^ is the sound velocity, the minimum 
value of the field f(z, t) is greater than zero and the wave is called gray soliton. If v = 
(stationary dark soliton) the minimum of the field f(z,t) is zero and the wave is called 
black soliton jTUEj- The velocity c s is the maximal velocity of the dark soliton, namely 
the velocity of a wave of infinitesimal amplitude, i.e. a sound wave, propagating in a 
medium of density = 0^. 



3. 3D bright and dark solitons 

We have seen that for the ID CNLSE the axial density profiles of bright and dark 
solitons do not depend on the sign of 5d and 5n, but only by their product 5d • o~n 
which discriminates between bright and dark solitons. It is important to stess that this 
property is valid only in the ID limit. In fact, as previously shown, the ID CNLSE is 
only the weak-coupling limit of ID NPSE, that is derived from the 3D CNLSE. In the 
remaining part of the paper we show that, in general, for a fixed value of the strength 
g, 3D bright solitons (Sp ■ Sn — 1) with anomalous dispersion (5d = 1) and self-focusing 
nonlinearty (Sn = 1) do not have the same axial density profile of bright solitons with 
normal dispersion (Sjj = —1) and self-defocusing nonlinearity (5n = —1). The same 
happens for 3D dark solitons. In the remaining part of the paper we call 3D solitons the 
solitary solutions of the NPSE, which are good approximations of the exact solutions of 
the 3D CNLSE under transverse confinement. 

To analyze 3D bright and dark solitons we start from the NPSE and set 

f(z, t) = (j)(z - V t)e i0{z - Vt) e i{v2/{2SD) -^ , (9) 

where ( = z — vt is the comoving coordinate of the soliton and both and 6(() are 
real fields. In this way we get an equation for 0(C), namely 

(* + J\ t = Jr. (f - e^) + , (io) 



25 D J 2 v > y/l — 5 N g(p 

and also an equation for the phase 0{(), that is 

vtf = 5 f (0 9" + 2<f>'e') . (11) 

Note that the prime means the derivative with respect to (; moreover, if = then 
from Eq. (11) one finds 6' = v/8r>. 

The two equations (10) and (11) are coupled but the second one can be written, 
for ^ 0, as 

M = iD iqi , (12) 

from which we find 

v 2 = 5 D (fO' + e , (13) 
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FIG. I. Axial density profile p(() = 0(C) 2 and transverse width a(() of the bright 
soliton (5 D • S N — 1). Strength of the nonlinearity: g = 1. A-F means 3D bright soliton 
with anomalous dispersion (5d = 1) and self-focusing nonlinearity (Sn = 1); N-D means 
3D bright soliton with normal dispersion (Sp = — 1) and self-defocusing nonlinearity 
(S N = -1). 



where the integration constant £ = for bright solitons and £ 
with 0oo the value of the field <f)(() at £ = ±oo. Because C 
Eq. (13) can be expressed as 



1 roo 



0^ for dark solitons, 
for bright solitons, the 



(14) 



By using this formula, the equation of 0(C) becomes 



^ =~90 ' 



where 



V(0) 



1 



v 2 4 



/i (j) 2 - 0Vi - 5jv g<P 2 + 



'D 



2 



(15) 



(16) 



Thus, the field 0(C) can be thought as the coordinate of a fictitious particle at time 
C- In this picture V(<j>) is the external potential acting on the fictitious particle. The 
constant of motion is then given by 



K 



\<P' 2 + v(<p) 



and after separation of variables we get 

rtrh 

'2 (K-V (</>)) 



(17) 



(18) 
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We stress that in our previous papers [13b] and ^3] only anomalous dispersion (5d = 1) 
and v (f 2 ^ = were considered. Here we investigate for the first time both normal 
dispersion (5d = — 1) and v (p 2 ^ ^ 0. In particular, the Eq. (14) and the centrifugal 
term of Eq. (16) are new and essential to study the case v </>^ ^ 0, i.e gray solitons with 
normal or anomalous dispersion. 

3.1. Bright solitons 

The NPSE admits bright solitons for 5d ■ 6~n = 1. In this case we have (j)^ = and 
K = 0. By using the Eq. (18) we find that the bright soliton <f)(() satisfies the implicit 
formula 



where ArcTanh[x] and ArcTan[x] are the hyperbolic arctangent and the arctangent, 
respectively. They satisfy the formula ArcTan[ix] = i ArcTanh[x] which is useful when 
Sd — — 1. It is not difficult to show that this equation reduces to Eq. (7) in the ID 
weak-coupling limit g<p 2 <C 1 taking into account the relationship between /x and g. This 
relationship is obtained by imposing the normalization condition to the bright soliton. 
Choosing / (p(() 2 d( = 1 we obtain 



It is clear that for 5d = — 1 it must be \x > 1, while for 5d = 1 one has < /x < 1 and 
the 3D bright soliton is dynamically stable only for 1/2 < /x < 1 [T3*| I14j. 

In Fig. 1 we plot the bright soliton with nonlinearity g = 1. The figure shows that 
the ID bright soliton given by Eq. (7) and the two 3D bright solitons given by Eq. (18) 
have different axial profiles because their transverse profiles are quite different. For the 
ID bright soliton a(() = 1, while for the two 3D bright solitons a(() is not constant and 
depends on the sign 5n of the nonlinearity: o~(Q > 1 with 5n = — 1 (self-defocusing and 
normal dispersion) and a(Q < 1 with Sn = 1 (self-focusing and anomalous dispersion). 

In Fig. 2 we plot the axial width Zp and transverse width cr(0) of the bright soliton 
as a function of the strength g. As expected, the axial width decreases by increasing 
g but for the A-F bright soliton, i.e. the 3D bright soliton with anomalous dispersion 
(6d = 1) and self-focusing (5n = 1), there is a critical strength g c = 4/3 = 1.33 above 
which the soliton does not exist anymore (see also 0). This critical value is in excellent 
agreeent with the numerical result g c = 1.35 obtained by solving the 3D CNLSE [21 J . 

For the N-D bright soliton, i.e. the 3D bright soliton with normal dispersion 
(8d = —1) and self-defocusing (5n = —1), there is not a critical strength and the 
N-D bright soliton exist for any g, as the strictly ID bright soliton. 





(19) 




(20) 
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FIG. 2. Axial full width half maximum zf and transverse width a at ( = as a function 
of the strength g of the nonlinearity for the bright soliton (5 D • S N — 1). The meaning 
of labels is the same of Fig. 1. 



The lower panel of Fig. 2 shows that the transverse width of the N-D bright soliton 
grows indefinitely and a remarkable consequence is that the soliton geometry changes 
from cigar-shaped to disk-shaped by increasing the strength g of the cubic nonlinearity: 
at g = 11.94 we find zf = o. 

3.2. Black and gray solitons 

The NPSE admits dark solitons for 5 D ■ 5n = —1- In this case we have (f)^ ^ and 
K = l^(0oo). By using the Eq. (17) we first numerically determine the dark soliton </)(() 
with v = 0. As previously stressed this dark soliton is called black soliton because its 
minimum value is zero. For the black soliton from Eq. (11) one finds 9' = 0. 

In Fig. 3 we plot the axial profile p(() = 0(C) 2 and the transverse width a(() of 
the black soliton setting g = 1 and 0oo = 1. We compare the ID black soliton given 
by Eq. (8) and the 3D black solitons obtained by numerically solving Eq. (15) with 
Eq. (16). The transverse width is a ~ 1 near the black hole but for 3D dark solitons 
it is cr = (1 — Sn^oo) 1 ^ & t infinity. It is clear that this condition bounds the domain 
of existence of the 3D dark soliton with self- focusing nonlinearity (5 N = 1): this dark 
soliton exist only for < 1. 

When v ^ then the dark soliton is called gray soliton and its minimum value is 
greater than zero. In Fig. 4 we plot the hole of the 3D gray soliton with anomalous 
dispersion (5f> = 1) and self-defocusing nonlinearity (5n = 1) for some values of v. By 
increasing the velocity v the depth of the hole is reduced and becomes zero when the 
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FIG. 3. Axial density profile p(() = </>(C) 2 and transverse width ct(£) of the dark 
soliton (Sd ■ $n — — 1) with v = (black soliton). Strength of the nonlinearty g = 1 
and asymptotic field = 1. A-D means 3D black soliton with anomalous dispersion 
(5d = 1) and self-defocusing nonlinearity (Sn = — 1); N-F means 3D black soliton with 
normal dispersion (5d = — 1) and self-focusing nonlinearity (5n = 1). 



velocity reaches the sound velocty c s , given in general by the formula 



\ 



$NgPc 



1 S N gp c 
4 n - 



\3/2 



(21) 



4 sfl - tfjv^poo 4 (l-^p c 

where = (p^ [22]- The 3D gray soliton with 5n = — 1 exist for any g if i> is 
sufficiently small. Nevertheless, it has been shown [22] that these 3D dark solitons 
become dynamically unstable for a large nonlinearity g (g = 1.5 for the black soliton) due 
to the so-called snake instability The snake instability implies the apparence of a non- 
axisymmetric purely imaginary eigenvalue in the Bogoliubov spectrum of elementary 
excitations. In this case the black soliton decays into a solitonic vortex [23]. 

In addition, it has been shown that many complex modes appear in the 3D black 
soliton at g — 4; for this nonlinearity an axisymmetric vortex ring emerges with energy 
intermediate between that of the black soliton and the solitonic vortex [23]- 



4. Conclusions 

We have studied 3D bright and dark solitons with anomalous and normal dispersion 
taking into account their transverse width by using a non-polynomial Schrodinger 
equation. We have found new general analytical solutions for 3D bright solitons under 
transverse harmonic confinement. In addition, we have also analyzed 3D black and 
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FIG. 4. Axial density profile p(() = 0(C) 2 of the gray soliton (v ^ 0) with anomalous 
dispersion (Sd = 1) and self-defocusing nonlinearity (5n = —1)- Strength of the 
nonlinearty g — 1 and asymptotic field 0^ = 1. 



gray solitons showing that self-focusing dark solitons (and with normal dispersion) 
exist only below a critical axial density. As a by-product we have obtained a formula 
for the sound velocity, i.e. the maximal velocity of propagation of a gray soliton, 
useful for both self-focusing and self-defocusing media. Atomic matter waves of Bose- 
Einstein condensates are now rutinely produced with dilute gases of alkali-metal atoms 
at ultralow temperatures. Also optical pulses in graded-index fibers are now available. 
As previously stressed, our predictions on 3D solitons under transverse confinement can 
be experimentally tested using these materials. 
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